Abstract. We present the X-ray spectrum of 
MR Vel /RX J0925.7— 4758, obtained with the Medium 
Energy Grating spectrometer of the Chandra X-ray 
Telescope. The simplest models used by earlier authors, 
stellar atmospheres in combination with a thermal plasma 
in collisional ionization equilibrium, cannot explain the 
spectrum. Neither does a photo-ionized plasma. We iden- 
tify P Cygni profiles of FeXVII and OVIII, from which 
we conclude that these lines arise in a wind. We conclude 
- - that major uncertainty exists about the bolometric 
. luminosity of MR Vel, and perhaps of supersoft sources 
' in general, so that the theoretical prediction that this 
£S) ■ luminosity derives from steady nuclear burning cannot be 
verified. 
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1. Introduction 

Supersoft X-ray sources are called such because they emit 
most of their X-ray flux at energies less than 0.4 keV. First 
discovered with the Einstein satellite (Long et al. 1981), 
they were recognized as a class only after the more sensi- 
tive observations with ROSAT (Trumper et al. 1991). The 
suggestion that the X-ray emission of supersoft sources is 
caused by steady nuclear burning at the surface of white 
dwarfs with masses near the Chandrasekhar limit (Van 
den Heuvel et al. 1992, see also Shara et al. 1977, Ibcn 
1982) gained credibility when the nova GQ Mus was found 
to be a supersoft source (Ogelman et al. 1993). 

Van den Heuvel et al. (1992) argued that the high 
mass-transfer rates (<; 10 _7 M Q yr _1 ) demanded by the- 
ory for steady-state hydrogen burning require donors with 
masses in the range of 1.4-2.2 M Q . However, optical iden- 
tifications soon showed that the class of supersoft X-ray 
sources is far from homogeneous, in some cases exclude 
k5 1.4 Mq donors, and also cast doubt on the premise 
that the white dwarfs in these systems accrete steadily at 
high rates (e.g. the review by Gansicke et al. 2000). 

Black body fits indicated that the X-ray luminosities 
of supersoft sources - several of which are located in the 
Magellanic Clouds, and thus have known distances - are 
in excess of the Eddington limit for a white dwarf (Grcincr 
et al. 1991). Fitting LTE white-dwarf atmosphere model 
spectra to the X-ray data, Heisc et al. (1994) showed that 
the X-ray luminosities may well be at or below the Ed- 
dington limit. NLTE models may give very different lu- 
minosities again (Hartmann & Heise 1997). The rather 
limited spectral resolution of the Einstein and ROSAT X- 
ray data prevents an accurate determination of the X-ray 
spectrum and of the interstellar absorption, and this leads 
to very large uncertainty of the estimated bolomctric flux. 
Accurate determination of the luminosities is important 
for comparison with the models, and to establish the re- 



lation of the supersoft sources with other sources showing 
very soft spectra at low luminosities, such as the variable 
source in the globular cluster NGC5272 whose nature is 
unknown (Dotani et al. 1999) or VSge whose X-rays may 
arise from colliding winds (e.g. Wood & Lockley 2000). 

Studies of Cal87 and MR Vel / RX J0925.7— 4758, 
both supersoft sources with significant count rates above 
0.4 keV, have shown that an LTE or NLTE white-dwarf 
atmosphere model does not correctly describe the X-ray 
spectrum (Hartmann et al. 1999, Ebisawa et al. 2001). 
In this paper we present the first X-ray spectrum of the 
supersoft source MR Vel with sufficient resolution to test 
spectral models in some detail. 
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Table 1. Summary of models for the X-ray data of MR Vel, giving the instrument used, its wavelength range and 
the resolution at IkeV, the model used and its fit parameters. BB, SA and TP indicate a black body, NLTE stellar 
atmosphere model, and thermal plasma in collisional ionization equilibrium, respectively. An index 2 indicates the 
second spectral component. Luminosities L are bolometric for an assumed distance of 1 kpc; column densities TVh are 
the same for both components, where applicable. References a. Motch et al. 1994 b. Hartmann et al. 1999 c. Ebisawa 
et al. 2001 d. this paper 





instrument 


E 


AE/E 


model 


T 


9 


L 


T 2 


U 




x 2 /dof 






(keV) 


(1 keV) 




(10 5 K) 


(cms -2 ) 


(ergs -1 ) 


(10 s K) 


(ergs -1 ) 


(cm" 2 ) 




a 


ROSAT PSPC 


0.1-2.4 


0.4 


BB 


5.6 




510 38 






1.7 10 22 


1.1/17 


b 


SAX LECS 


0.1-10. 


0.15 


SA 


8.75(5) 


10 9 


4.1 10 36 






1.5(2) 10 22 


5.6/15 


b 


SAX LECS 


0.1-10. 


0.15 


SA+TP 


8.75(4) 


10 9 


2.9 10 36 


80(60) 


< 1.5 10 35 


1.4(2) 10 22 


3.4/14 


c 


ASCA SIS 


0.4-10. 


0.08 


BB 


5.2 




1.8 10 40 






2.1 10 22 


8.1/73 


c 


ASCA SIS 


0.4-10. 


0.08 


SA+TP 


10.9 


10 10 


710 35 


12 




1.3 10 22 


2.0/71 


(1 


Chandra HETG 


0.4-8.0 


0.02 


BB 


5.2 




2.8 10 39 






1.6 10 22 


2.5/252 


d 


Chandra HETG 


0.4-8.0 


0.02 


SA 


8.97 


10 9 


5.3 10 35 






1.0 10 22 


2.9/252 


d 


Chandra HETG 


0.4-8.0 


0.02 


SA+TP 


8.96 


10 9 


4.1 10 35 


20.6 


4.9 10 34 


1.0 10 22 


2.6/250 



2. Earlier observations of MR Vel 

MR Vel was discovered in the Galactic Plane in the 
ROSAT All Sky Survey, in a search for supersoft sources. 
The intrinsic softness of the source combined with strong 
interstellar absorption cause most of the photons of this 
source to be detected between 0.5-1.5 keV. A black body 
fit to the data from a ROSAT PSPC pointing gives a tem- 
perature of ~ 48 eV, for an interstellar absorption column 
of - 1.7 10 22 cm -2 (Motch et al. 1994). Observations with 
the BeppoSAX MECS and the ASCA SIS show excess flux 
at energies <; 1.2 keV with respect to the flux predicted 
by black body or stellar atmosphere models. To explain 
this, a thermal plasma has been suggested as a second 
component of the spectrum. However, models combining 
a NLTE white-dwarf atmospheric model with a thermal 
plasma in collisional ionization equilibrium do not provide 
acceptable fits to the data, which appear to show absorp- 
tion edges at 0.87keV (OVIII), 1.02keV (unidentified), 
and 1.36 keV (NeX), according to Ebisawa et al. (2001). 

NLTE atmosphere models give a much lower (unab- 
sorbed) flux for the observations than LTE models, which 
in turn give lower fluxes than black body fits. Whereas 
LTE models can give near-Eddington luminosities for 
MR Vel at a distance of 1 kpc, NLTE models lead to much 
larger distances or lower luminosities. 

Spectral fits to pre-Chandra X-ray data are summa- 
rized in Table |l|. 

A VI th magnitude optical counterpart was identified 
on the basis of a ROSAT HRI position (Motch et al. 
1994) and showed a photometric and spectroscopic period 
of 4.029d (Schmidtke et al. 2001). The optical lightcurve 
is asymmetric and single-peaked, with a full amplitude 
of about 0.3 mag. A shorter, unexplained period of about 
0.25 d also appears to be present in the photometry. The 
very red optical continuum - presumably of the accretion 
disk - and the very strong interstellar absorption lines in- 
dicate a high interstellar absorption, in agreement with 



the X-ray results (Motch et al. 1994). Transient satellite 
lines of the Ha line indicate a transient double-sided jet, 
with a velocity of ~ 5000 km s -1 , and an opening angle 
between 34 and 80 degrees (Motch 1998). The photometry 
shows no eclipses, indicating an inclination i <, 65°. 

The radial velocity curve of the He II 4686 and Her 
lines, with a semi-amplitude of 75kms _1 , with the ab- 
sence of eclipses lead to a likely donor mass between 1.0 
and 1.6 M for assumed masses between 0.5 and 1.4M 
for the white dwarf. The absence of the donor in the visual 
spectrum indicates a bright disk, which in turn indicates 
a minimum distance of a kpc. 

3. Chandra observations and data analysis 

MR Vel was observed on 14/15 November 2000, for a net 
exposure time of 56.5 ks, using the High-Energy Trans- 
mission Grating Spectrometer (HETGS) of the Chandra 
X-ray Observatory (Canizares et al. 2000 and references 
therein). The high energy grating (HEG) could not be 
used for analysis because of a low S/N-ratio. The medium 
energy grating (MEG) is optimized for the range of 1.2- 
30 A and provides a spectral resolution of 0.02 A FWHM. 
The photons are collected at an array of CCDs, ACIS-S, 
in the focal plane. 

The data were processed with standard Chandra X- 
ray Center software CIAO-2.1, as follows. Photons are ex- 
tracted by applying a spatial filter around the spectral 
image, and negatively and positively dispersed spectra are 
added. Most background photons are recognized by pulse 
height analysis, and removed; the higher-order spectra of 
this faint source do not contain useful information. Flux 
calibration was done with the 1999-07-22 version of the 
detector response matrix, part of the Chandra calibration 
database CALDB 2.3. 

The errors in each bin with N photons are approxi- 
mately Poissonian and given by (see Gehrels 1986): 
AJV = 1 + y/N + 0.75 (1) 
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ion 


number of 


quantum 


number of 




levels 


number 


transitions 


HI 


5 


5 


- 


Hell 


3 


3 


- 


CV 


17 


4 


315 


CVI 


10 


4 


285 


NVI 


17 


4 


315 


NVII 


10 


4 


285 


OVI 


5 


3 


153 


VII 


17 


4 


315 


VIII 


10 


4 


285 


NeVII 


6 


2 


2111 


NeVIII 


5 


3 


153 


Ne IX 


17 


4 


315 


NeX 


10 


4 


285 


FeXV 


19 


3 


5074 


FeXVI 


7 


4 


184 


FeXVII 


21 


3 


3409 


Fe XVIII 


2 


2 


9437 


FeXIX 


6 


2 


23915 


FeXX 


8 


2 


37795 


FeXXI 


12 


2 


34017 



Table 2. The ions used for the construction of the model 
atmospheres in TLUSTY, with the respective number of 
energy levels and the corresponding quantum number up 
to which all levels are considered. The number of line tran- 
sitions used by SYNSPEC are also given; these are all 
transitions between energy levels with quantum numbers 
from 1 to 10. 



Models for thermal plasmas in collisional ionization 
equilibrium were computed in SPEX, with the code de- 
veloped and described by Kaastra et al. (1996). 



For comparison with model spectra, we use x 2 fitting and 
bin the data, until the oversampling of the spectral reso- 
lution is less than a factor 3, and until at least 20 counts 
are present in each energy bin. The resulting spectrum is 
shown in Fig. [H 



4. Model spectra 

We compute NLTE model spectra for white dwarf at- 
mospheres, characterized by effective temperature, grav- 
itational acceleration, and elemental abundances, as de- 
scribed in Hartmann et al. (1999), i.e. using the codes 
TLUSTY (version 195) for the construction of a model 
atmosphere, and SYNSPEC (version 42) for the compu- 
tation of the emergent spectrum. These codes were devel- 
oped and are described by Hubeny (1988) and Hubeny & 
Lanz (1995). For the atmosphere model we use 70 depth 
zones, and a frequency range from 10 12 Hz to 9.610 17 Hz; 
convergence is considered to have been achieved when the 
maximum relative change between subsequent model pa- 
rameters is less than 0.001. The ions and transitions used 
in our computations are listed in Table]^; their numbers 
are limited by computational expense. The atomic data 
are extracted from the Opacity Project with TOPbase (see 
Cunto et al. 1993). 
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Fig. 1. The observed X-ray spectrum of MR Vel binned to the minimum width required for either 1/3 of the instru- 
mental spectral resolution or 20 counts, together with the best black body fit (dotted line) and the best fit combination 
(thick line) of the spectra from a NLTE atmosphere and a thermal plasma in collisional ionization equilibrium (the 
latter is also shown separately with a thin line). Some identified spectral features are indicated. 



5. Results 

The observed spectrum shown in Fig.[l] shows the fol- 
lowing emission lines (with theoretical wavelengths from 
Phillips et al. 1999). The lines at A 15.014, 15.265, 16.780, 
17.055, and 17.100A are all due to FeXVII. The line at 
A 18. 973 A is the OVIII Lyman a line. The strong line 
complex between A 14.3-14.6 A is a blend of Fe XVIII lines, 
which cannot be identified individually, the strong feature 



at A 14.6-14.7 A may be a blend of O VIII and Fe XIX lines. 
The complex of the NelX A 13.448, 13.553 and 13.700A 
triplet, respectively the resonance line (w), the intercom- 
bination (x+y) lines, and the forbidden (z) line of the 
He-like ion; these lines are not recognizable individually, 
but their presence is clear from the enhanced flux in their 
wavelength range. 

We have fitted the observed spectrum with the models 
discussed by previous authors, i.e. black body, LTE atmo- 



Henk Bearda et al.: X-ray spectrum of MR Vel 



6 

sphere and NLTE atmosphere spectra. Following earlier 
authors, we also added a thermal plasma in collisional ion- 
ization equilibrium to a NLTE atmosphere. In Fig. |l| we 
show the spectra of the best fit black body and the com- 
bination of the atmosphere and the plasma. The model 
parameters are given in Table |l]. 

The interstellar absorption edges due to Nel (at 
A 14.30 A) and Fe-L (at A 17.54 A) are recognizable in the 
models and in the data, and confirm that interstellar ab- 
sorption is appreciable. It is clear that neither model is ac- 
ceptable. The observed spectrum shows much more struc- 
ture than a black body; but not the structure expected 
for a calculated atmospheric spectrum. This is true for a 
gravitational acceleration of 10 9 cms -2 and effective tem- 
peratures between 8.00-9.15 10 5 K. 

The combination of a thermal plasma with the NLTE 
atmosphere indeed adds lines, but the strengths of these 
are well below those of the observed emission complexes 
between A 14-16 A. The flux contribution of the plasma is 
so small that the parameters of the atmospheric model 
are virtually the same for both models with and without 
an added thermal plasma. In some calculations we allow 
an arbitrary global shift in the wavelengths of the lines 
from the thermal plasma to take into account the current 
accuracy of the HETG-MEG absolute wavelength calibra- 
tion; this doesn't change any of the conclusions. Similarly, 
varying the abundances of iron or oxygen in the thermal 
plasma only results in marginal improvement; only some 
emission lines can be reproduced better this way. Assum- 
ing different interstellar absorption densities for the two 
model components, complicates the analysis further with 
no better results. 

If one tries to enhance the emission line strengths 
by changing the temperature or the emission measure of 
the plasma, one always finds very strong lines at wave- 
lengths where none are observed. For example, whenever 
the strength of the FcXVII lines at A 15.014 and 15.265 A 
are enhanced to bring their fluxes to the observed level, the 
lines from the same ion at A 12.124 and 12.264 A become 
too strong. It also produces a thermal continuum which 
is significantly higher between A 7-12 A than the observed 
flux. 

To compare the relative strengths of the identified 
FeXVII lines with theoretical predictions, we have fitted 
them with Gaussian profiles, as printed in Table||. No 
thermal plasma in collisonal ionization equilibrium can 
reproduce the observed relative strengths of the FeXVII 
lines, even when we allow a continuous distribution of 
emission measures as a function of the temperature (the 
differential emission measure model, see Liedahl et al. 
2001). Also the strength of the Fe XVIII lines between 
A 14. 3-14. 6 A can't be explained by such plasma. These 
results remain true when the abundances are varied from 
solar, because they affect the relative strenghts of lines 
from different elements, not those from the same element. 



Theoretical 
wavelength (A) 


Fitted 

wavelength (A) 


Relative flux 


15.014 
15.265 
16.780 

17.055/17.100 


15.05 ±o.oi 
15.276 ±o.oo7 
16.79 ±o.oi 
17.10 ±o.oi 


0.63 -0.22 / +0.25 
1.00 -0.24 / +0.27 
1.61 -0.7S / +1.40 
3.45 -1.44 / +2.20 



Table 3. Parameter values for the fit of the identified 
FeXVII lines with gaussian profiles. The two lines at 
A 17.055 and 17.100 A are fitted with one profile. The 
fluxes are printed relative to the flux of the line at 
A 15.265 A, and have been corrected for interstellar ab- 
sorption with an absorption measure of 0.99 1 22 cm~ 2 . 
Note that this emission feature may have a contribution 
of the O VIII Lyman 7 line. 



The main problem is the ratio of the 15.014 / 15.265 
lines, for which optically thin models predict a value of 
3-4 for any density or temperature, while we observe 
0.63. The ratio of the 16.780 and 17.055/17.100 lines 
to the 15.265 line is however consistent with an optically 
thin plasma. Dropping the optically thin assumption, the 
15.014 line can be brought into agreement with the 16.780 
and 17.055 / 17.100 fluxes for pure FeXVII column densi- 
ties of the order of 10 19 cm -2 , but also in that case the 
predicted 15.014/ 15.265 ratio still remains similar as in 
the optically thin case. Thus, also for the optically thick 
model the 15.265 line is 6 times weaker than observed. 

We conclude that none of the earlier models describe 
the observed X-ray spectrum of MR Vel. 

We investigate whether the spectrum shows signa- 
tures of a photo-ionized plasma, as listed by Liedahl et 
al. (2001), viz. a) narrow radiative recombination con- 
tinua; b) large ratio (x + y + z)/w of the He-like triplet; 
c) the weakness of L-shell iron emission compared to K- 
shell emission from lighter elements. The strength of the 
Fe XVII lines and the relative weakness of the Ne IX lines 
argues against a photo-ionized plasma. The observed spec- 
trum also does not show evidence for narrow recombina- 
tion continua. The strengths of the individual Ne lines in 
the spectrum are too small to determine their ratios. We 
conclude that a simple photo-ionized plasma also doesn't 
describe the observations. 

6. P Cygni profiles indicate a stellar wind 

Considering that none of the previously suggested mod- 
els describes the data, we have another look at the ob- 
served spectrum. Fig. ^| gives an expanded view of the five 
strong, individually identifiable lines and the spectrum 
around the OVIII Lyman j3 line at 16.003 A. The simi- 
larity between the profiles of these lines indicates that the 
P Cygni structure seen in them is real. One would expect 
to see a similar structure for the Lyman 7 and S lines; in- 
spection of the spectra near their wavelengths (15.175 and 
14.817 A, respectively) indeed shows indication of redward 
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emission and blueward absorption. We conclude from the 
P Cygni profiles that the stronger individually identifiable 
lines (see Sect. 5) in our spectrum arise in a wind. 

The absorption parts of the OVIII Lyman a and (3 
and FeXVII A 17.055 A lines are saturated; the absorption 
part of the O VIII Lyman S line is not saturated. For the 
remaining lines the situation is ambiguous. 

Strong P Cygni lines are present in the spectra of high- 
mass stars, e.g. Wolf-Rayet stars, and also in the spectra of 
systems with an accretion disk. The radial velocity curve 
and the absence of strong stellar lines in the optical spec- 
trum, exclude the presence of a high-mass star in MR Vel 
(see Sect. 2). Thus, the P Cygni profiles in the X-ray spec- 
trum of MR Vel are more likely related to the presence of 
an accretion disk. 

In cataclysmic variables, the terminal wind velocities 
are often fairly high, reaching values of 5000 km s -1 (e.g. 
Woods et al. 1992). The similar velocity of the transient jet 
seen in the optical has been used by Motch (1998) to argue 
that the accreting object is a white dwarf. Remarkably, 
the highest absorption velocities seen in Fig. || are rather 
lower, of the order of 1500 km s -1 . 

It seems possible that a combination of an atmosphere 
spectrum with P Cygni profiles from a wind gives a good 
description of the X-ray spectrum of MR Vel. It is beyond 
the scope of this article to construct such a model. At 
this point we repeat our earlier statement that no obvious 
atmospheric features can be seen in the X-ray spectrum, 
as shown in Fig. |l|. 

7. Discussion 

The observed X-ray spectrum of MR Vel is dominated by 
emission features, and is not the continuum with absorp- 
tion features expected on the basis of stellar atmosphere 
models, nor is it an optically thin spectrum of a thermal 
plasma in collisional ionization equilibrium, or a photo- 
ionized plasma. Instead, we suggest that the spectrum is 
a combination of a stellar atmosphere and a wind. 

The bulk of the luminosity of the models fitted by ear- 
lier authors lies outside the observed X-ray regimes, at 
softer energies. As a result the luminosities derived from 
the same observations range from above the Eddington 
limit for a black body to less than 0.001 of the Eddington 
limit for the spectrum of a NLTE atmosphere, as can be 
observed in Table [I]. The uncertain magnitude of the inter- 
stellar absorption adds to the uncertainty. This also may 
be seen in Table [l], by comparing the same model fitted to 
data of different satellites. 

Therefore we cannot predict the bolometric luminos- 
ity corresponding to the X-ray spectrum of MR Vel. This 
remains true for models with a wind. We can only de- 
termine the luminosity in the observed range of 8-20 A 
at ~ 10 35 ergs _1 for an assumed distance of lkpc. As 
a result, we can no longer be certain that MR Vel is as 




Fig. 2. Profiles of some strong, individually identifiable 
lines in the spectrum of MR Vel. The average flux in each 
frame has been normalized to unity; this level is indicated 
with the dashed lines. The wavelength scale is given in 
Doppler velocities, based on central wavelengths shown on 
the right of each frame; the wavelength resolution is the 
instrumental resolution, i.e. more points are shown than 
in Fig.[j]. The central wavelengths of the oxygen Lyman 
7 and 5 lines are indicated with vertically dotted lines in 
the two lower frames. 
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luminous as expected for a white dwarf steadily burning 
helium at its surface. 

One of the other supersoft sources, Cal 83, has been ob- 
served with the reflection grating spectrograph of XMM, 
with a resolution of 0.05 A in the range of interest for 
this softer source, 20-38 A (Paerels et al. 2001). At first 
glance, the observed spectrum of Cal 83 could consist of 
either emission features with a width of ~ 1 A, or a contin- 
uum with many absorption features of comparable width. 
Paerels et al. prefer the latter interpretation, but cannot 
match the absorption features with features from an at- 
mosphere dominated by CNO, or from a model with ele- 
ment abundances as appropriate for the LMC, up to Fe. 
At the moment then, the bolometric flux of Cal 83 also is 
unknown. 

Before definite conclusions can be drawn about the na- 
ture of the supersoft sources, more progress is required in 
the theoretical description of their X-ray spectra. In par- 
ticular, it seems necessary to develop models that combine 
atmosphere models with winds. The atmosphere models 
themselves can be improved with respect to the current 
models in several ways. More elements and ions can be in- 
cluded in the computation of the atmospheric structure, 
with a more complete account of the line-blanketing ef- 
fects. Gradients in temperature and gravitational acceler- 
ation, as can be expected for accretion disk spectra, cannot 
be ruled out. Non-static or outflowing atmospheres may 
have to be considered. The question arises whether the 
jet observed in MR Vel and the winds observed in several 
supersoft sources contribute to their X-ray spectra. 
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